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Abstract 

We use the analytical model recently introduced in Ref. Q, to 
investigate the statistics of temperature fluctuations on the cosmic 
microwave background (CMB), induced by topological defects. The 
cases of cosmic strings and textures are studied. We derive analyti- 
cally the characteristic function of the probability distribution for ^ 
and use it to obtain the lowest twelve moments including the skewness 
and the kurtosis. The distribution function is also obtained and it is 
compared with the Gaussian distribution thus identifying long non- 
Gaussian tails. We show that for both cosmic strings and textures all 
odd moments (including skewness) vanish while the relative deviation 
from the Gaussian for even moments increases with the order of the 
moment. The non-Gaussian signatures of textures, derived from the 
distribution function and the moments, are found to be much more 
prominent than the corresponding signatures for strings. We discuss 
the physical origin of this result. 



1 Introduction 



Theoretical models for large scale structure formation can be divided in two 
classes according to the type of primordial perturbations they consider; mod- 
els based on adiabatic Gaussian perturbations produced during inflation and 
models based on non-Gaussian perturbations which are naturally provided 
by topological defects. The cold dark matter (CDM) model, based on Gaus- 
sian primordial fluctuations, produces an evolved density field with small and 
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intermediate scale structure in reasonable agreement with observations!^. 
Recent observations on large scales however, have created significant chal- 
lenges for the CDM model. One of the most serious such challenges comes 
from the recent detection of anisotropy in the CMB by the DMR (Differ- 
ential Microwave Radiometer) instrument of COBE (Cosmic Background 
Explorer) satellite 

This discovery has provided a new powerful experimental probe for test- 
ing theoretical models for large scale structure formation. The temperature 
sky maps constructed by DMR were used to obtain the rms sky variation 
^< {^f > (where AT = T{9i) - T{92), and 61-62 = 60° is the beam 
separation in the COBE experiment) and the rms quadrupole amplitude. A 
fit of the data to spherical harmonic expansion has also provided the angu- 
lar temperature auto-correlation function C{A6) =< ^{6)^{6') > where 
<> denotes averaging over all directions in the sky, 5T{6) = T(6)— < T > 
and A6 = 6 — 6'. This result was then used to obtain the rms-quadrupole- 
normalized amplitude Qrms-PS and the index n of the power law fluctuation 
spectrum assumed to be of the form P{k) ~ k"'. The published results are: 

n = 1.1 ±0.5 

Qrms-PS = (5.96 ± 0.75) X 10-*^ (1.1) 
AT 

i^)rms = (1.1 ±0.2) X 10-5 

Severe constraints are imposed on several cosmological models due to these 
results. For example, the CDM model with bias 1.5 < < 2.5 is incon- 
sistent with the COBE results for Hq > 50km/ (sec ■ Mpc) and is barely 
consistent for Hq ~ 50km/{sec • Mpc) Q It is therefore interesting 
to investigate the consistency of alternative models with respect to the 
COBE measurments. The natural alternative to models based on adia- 
batic Gaussian perturbations generated during inflation are models where 
the primordial perturbations are created by topological defects like cosmic 
strings global monopolesQ or textures Q. 

In a recent paper Q we introduced an analytical model and used it to 
study the effects of cosmic strings on the microwave background. Our model 
was based on counting random multiple impulses, inflicted on photon trajec- 
tories by the string network between the time of recombination and today. 
After constructing the temperature auto-correlation function, we used it to 
obtain the effective power spectrum index n, the rms-quadrupole-normalized 
amplitude Qrms-PS and the rms temperature variation smoothed on small 
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angular scales. For the values of the scaling solution parameters obtained 
in Refs.[^],[^ we showed that 



n = 1.14 ±0.5 

Qrms^PS = (4.5 ± 1.5)G// (1.2) 
AT 

( rp )rms — 5.5G'/i 

where /i is the mass per unit length of the string (the single free parameter in 
the cosmic string model for structure formation) and G is Newton's constant. 
Demanding consistency of our results with the COBE results (1.1) leads 

toi 

= (1.7±0.7) X 10-*^ (1.3) 
in good agreement with direct normalization of n from galaxy ]lO| and large 



scale structure [11] 1 12 1 observations (for more recent studies of the cosmic 
string model for structure formation see Ref. fl^). We concluded that 
the cosmic string model remains a viable model consistent with the up to 
now announced COBE data. Similar results to those presented in Ref. 
have been obtained by numerical simulations of the string network from 



the time of recombination to a small redshift|14|. These studies however, 
in contrast to our analytical model, do not attempt to take into account 
compensation (for a recent study of the effects of compensation see Ref. 
]l5[) and allow the string deficit angle to extend over the whole volume of 
their simulation. In addition they are constrained to fixed values of the 
scaling solution parameters produced by simulations which are still subject 
to some controversy [|^] [|^] § . Other interesting studies |0 have used the 
old picture for the cosmic string network]^], based on loops, to analytically 
calculate the effects of strings on the CMB. Recent simulations!^ @ however, 
have shown that the dominant component of the scaling solution consists of 
long strings rather than loops. 

Our analysis had utilized the data concerning the amplitude and the 
spectrum of the detected fluctuations in order to test the cosmic string model. 
This type of test can check the consistency of the model with the data but 
it can not distinguish it from other consistent theories. There are two basic 
reasons for this: First, the n = 1 Zeldovich spectrum is fairly generic in 
physically motivated theories and second, other theories like standard CDM, 
can also pass the amplitude normalization test ||l^ by utilizing tensor mode 
perturbations. 
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It is therefore clear that further tests are needed in order to distinguish 
between the topological defect models and other theories. One of the most 
interesting such tests is the subject of this paper: the study of the statistics 
of the CMB fluctuations. This is a particularly interesting issue in view of 
the existing temperature fluctuation sky maps obtained by COBE which 
are currently subjected to statistical analysis by the COBE collaboration. 
Such an analysis can reveal the characteristic signature of the source that 
produced the CMB fluctuations. The identification of this signature for 
topological defect models is the focus of the present work. 

Models based on fluctuations generated during inflation predict (in their 
generic form) a Gaussian distribution function for the CMB temperature 
fluctuations. On the other hand, models based on topological defects, like 
cosmic strings or textures, are distinguished due to the particular type of 
non-Gaussian fluctuations they predict. 

There is an extensive literature on the statistical effects of various types 
of non-Gaussian perturbations (for a small sample see Ref. |2^). On the 
other hand the literature on the statistics of seed-like perturbations is much 



more limited (see e.g. |^|Qp3|). Previous analytical studies on the type 
of non-Gaussian statistics induced by cosmic strings and other seed-like per- 
turbations have focused on the effects of point-like seeds on the statistics of 
large scale structure |2^. They have considered the superposition of over- 
dense, spherically symmetric kernels, thus obtaining the density distribu- 
tion function and other statistical properties. Those studies, even though 
elegant and convienient for large scale structure considerations can not be 
easily applied to the CMB case. The basic reason for this, is the need for su- 
perposition of variable size (and in general non-spherical) kernels to account 
for the effects of compensation in different Hubble times and horizon scales. 
The multiple impulse approximation of Ref. [l[ (see also Ref. |25| for an 
alternative application of the method) can incorporate this kernel variability 
and it is the method that we used in this work. 

In particular, we use the multiple impulse approximation to derive the 
characteristic function, the temperature distribution function and several 
moments of cosmic string and texture induced CMB fluctuations. These re- 
sults are then compared with the corresponding Gaussian results in an effort 
to find signatures of the cosmic string and texture models. In section 2 we 
give a brief review of our model in order to clarify the basic assumptions 
made (for a more complete account see Ref.[|l|). We also associate with it 
a set of statistical experiments. In section 3 we introduce some basic sta- 
tistical quantities and derive the characteristic functions for the statistical 
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processes described in section 2. In particular we derive the characteristic 
functions that correspond to string, texture and Gaussian temperature pa- 
terns. Finally, in section 4 we use the characteristic functions to obtain the 
probability distributions and the lowest twelve moments in each of the three 
cases. We also compare and discuss the obtained results. 



2 The Model 

We begin by reviewing the model introduced in Ref. which was originally 
designed to approximate the cosmic string induced temperature fluctuations. 
It will be seen that with some modifications it can also describe texture 
induced temperature fluctuations. 



According to cosmic string simulations 1 16 at any time t there are 



about 10 long string segments with a typical curvature radius t passing 



through each horizon volume. There is a well deflned mechanism[26| by 
which these segments give rise to localized linear temperature discontinuities 
(Fig. 1). Photons passing on different sides of a long straight string moving 
with velocity Vg perpendicular to the line of sight, reach the observer with 
a Dobbler shift pi] 

— = ±47rG^Us7s 

Our goal is to find the combined effects of temperature fluctuations induced 
by strings present between the time of recombination tj-ec and today. 
We approximate the photon path from the recombination time t^ec 

to 

the present time to by a discrete set of M = 16 Hubble time-steps ti such 
that tj+i = 2ti ~ 2^'^ for Zrec = 1400). We assume that the effects of 
long strings between the time of recombination and today give the dominant 
contribution to temperature fluctuations and therefore we consider photons 
emerging from the last scaterring surface at tree with a fixed uniform temper- 
ature. A beam of photons coming from a fixed direction will initially suffer 
rih 'kicks' from the Uh — 10 long strings within the horizon at ti = tree 
(the linear superposition of 'kicks' is consistent with the multi-string metric 
presented in Ref. [p7|| ) . Each temperature 'kick' from a string with arbitrary 
orientation with respect to the observer induces a fluctuation ^ of the form 

^ = ±47rG/x/3 (2.1) 

with 

(3 = k- {vs'js X s) 
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where k, s and Vg are the unit wave-vector, the unit vector along the string 
and the string velocity vector respectively. The sign changes along the 
string||2^. At the next Hubble time-step t2, Uh further 'kicks', uncorrelated 
with the initial ones, will be inflicted on the photon beam by the strings 
present within the new horizon scale t2 = 2ti. The process will continue 
until the M = 16 Hubble time-step corresponding to the present time to 
(Fig. 2). 

The physical process described above, corresponds to the following set 
of successive statistical experiments. For simplicity we will first consider the 
one dimensional case. The generalization to the realistic two dimensional 
case will then be made in a straightforward way. 

Consider a one dimensional (continuous or not) set of temperature pixels 
with initially uniform temperature distribution and fixed total length L. 
Consider now a localized step-function perturbation imposed on this surface. 
Such a perturbation would increase the temperature of a fixed small length Iq 
by an amount ^ = 6 but would decrease the temperature of a neighbouring 
equal length by the same amount (see Fig. 3). Thus after this 'trial' each 
pixel of the set has probability p(l) = ^ = pQ to have been positively 
perturbed, the same probability to have been negatively perturbed (p(2) = 
p(l) = po) and probability p{3) = 1 — ^ = 1 — 2po to have remained 
unperturbed. Let this 'trial' repeat no times before the first 'experiment' is 
completed. The next step is to repeat this experiment with a new scale for 
the step-function li = 21q and rii = no/2 number of 'trials'. The successive 
experiments continue until Im = ^o2*^ and um = no/ 2^'^. We demand 

= ■§ and riM = where n/j is a fixed positive integer to be identified 
with the number of long strings per horizon volume (see below). Therefore 
^0 = 2Aj+i and no = nh2^^ . This implies that 

Pj = 2^P0 = ^ (2.2) 

= ^ = nn— (2.3) 

with j = 0, 1, ...,M. The above described statistical process corresponds to 
our physical model provided that the following identifications are made: 

• The fixed length L is identified with the scale of the present horizon. 

• Each step-function perturbation with scale Ij is identified with a cosmic 
string perturbation induced at the Hubble time-step tj. At this time 
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the horizon scale is tj = 21 j. Compensation confines the string induced 
perturbation within the horizon scale 21 j. Clearly 21q is to be identified 
with the horizon scale at recombination while tM = 2Za/ (M is the final 
step) should be identified with the present horizon scale L. 

• Each 'experiment' is identified with a Hubble time-step. During the 
jth. 'experiment' there are nj = "''J 'trials' (string perturbations) 
corresponding to the number of strings per horizon length (volume) n/j 
times the total number of horizons within the fixed length L at the 
jth step. Clearly for j = M there is only the present horizon included 
and the total number of 'trials' (string perturbations) is njv/ = n^. 

The above described process needs to be improved in two ways in order to 
correctly approximate the physics: 

1. The one dimensional set of pixels must be promoted to a two dimen- 
sional one. 

2. The amplitude of the step-function perturbation must be allowed to 
vary in order to account for varying string velocities (varying param- 
eter f3 in equation (2.1)). Without introducing this improvement our 
analysis would still be interesting but it would approximate the pertur- 
bations induced by textures rather than strings. In the case of textures 
there is no velocity parameter involved but there are still equal mag- 
nitude positive and negative fluctuations depending on whether the 
photon 'falls in' or 'climbs out of the texture p^ . 

Before we incorporate these modifications we proceed to study the statistics 
of the perturbations in the above described simplest case. This will help 
illustrate our method more clearly while making its generalization a simple 
and straightforward task. 

3 The Statistics 

Let us first focus on a particular 'experiment' j (Hubble time-step). Define 
f(nj, kj) to be the probability that any given pixel will have been perturbed 
by ^ = kj6 at the end of the rij 'trials' (string perturbations). Since at 
each 'trial' there are three possible outcomes with known fixed probabil- 
ities p(l) = p{2) = = Pj, p{3) = 1 - 2^ = 1 - 2pj, the probability 
distribution /{rij, kj) may be obtained from the well known trinomial distri- 
bution (a simple generalization of the binomial). The trinomial distribution 
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fn-{x{l), x(2), x(3)) gives the probability that after rij 'trials' with three pos- 
sible outcomes, there are x{i) occurences of outcome i {i = 1, 2, 3) (obviously 
x(l) + x(2) + x{3) = rij). The trinomial distribution is: 

UM^),x{2),x{3)) = ^(i),^'^^j,^(3), P(l)"^S(2)-(^V(3)-(^) (3.1) 

Let x{l) (x(2)) be the number of positive (negative) temperature shifts for 
a given pixel while x{3) is the number of 'trials' that lead to no shift. Using 
the relations x(l) + a;(2) + a;(3) = rij and kj = x(l) — x(2) to change variables 
from rc(l), a;(2), x(3) to nj,kj,x{3) and summing over the possible x{3) we 
obtain f{nj, kj): 

_ nj\p]'-''^^\l - 2pjr^^) 

fin,, k,) - ^1,^ j^^TEEm^^^^EhEM^. ^^-'^ 

where the sum extends over all integer values of x{3) for which rij ± kj — x(3) 
is even, (notice that terms involving factorials of negative numbers vanish 
automatically) . 

A more useful and much simpler function that describes the statistics is 
the characteristic function (p{n, to) of the distribution. This has two impor- 
tant properties: 

1. It is the Fourier transform of the probability distribution i.e. 

n 

0(n,a;)= ^ e'^''f{n,k) (3.3) 

k=—n 

/(n,fc) = ^^'%— V(n,u;) (3.4) 

2. It can generate all moments < fc"* > of the distribution by differenti- 
ation i.e. 

< fc- >= (-l)™i-^0(„,a;)U=o (3.5) 

Using the property 

{p{l) +p{2) +pi3)r = E fn{xil),x{2),xi3)) 

x{i)+x{2)+x{3)=n 
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and equation (3.3) it is straightforward to show that the characteristic func- 
tion for the variable kj = x(l) — x{2) is: 



cj){nj,p,,u;) = {2pj cos(a;) + (1 - 2pj))"^ (3.6) 

However, we are interested in a multiple 'experiment' process i.e. we are 
looking for the distribution function of the variable: 

M~16 

K= (3.7) 

i=0 

where K is to be identified with the total temperature fluctuation ^ at 
the present time to- It is straightforward to show [^] that the characteristic 
function for a sum of independent random variables is equal to the product of 
the individual characteristic functions. Therefore, the characteristic function 
$(w) corresponding to the variable K is: 

M~16 
j=0 

This result may now be used to obtain the probability distribution by the 
Fourier transform (3.4). It may also be used to obtain all the moments of 
the distribution either by differentiation (using equation (3.5)) or by direct 
integration, using the distribution function. However, we must first identify 
correctly the parameters nj, pj for the physical problem under considera- 
tion. Clearly, equations (2.2), (2.3) need to be modified to account for the 
propagation of a surface (photon wavefront) rather than a line, through the 
string network. It is straightforward to generalize arguments leading to (2.2) 
and (2.3) to the two dimensional case to obtain: 

= 4Vo = 1^ (3.9) 



^0 r.im 



Using (3.9) and (3.10) in (3.8) we obtain: 

M~16 

$,(a;) = Yl (4^-^^ cos(a;) + (1 - 4i-A^))«'.4(«-^) ^^ -^-^-^ 

j=0 
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where the subscript t denotes that this result is vahd only for the case of 
texture-like perturbations where the magnitude of the individual pertur- 
bations (step-function magnitude) can be considered fixed. In the case of 
strings wc must generalize (3.11) further in order to account for the variable 
parameter (3 of equation (2.1) which represents the velocity and string ori- 
entation dependence of the perturbations. This generalization will lead to 
the multinomial distribution. 

Consider the above described process of successive perturbations with 
the additional feature of allowing Q possible magnitudes for the applied step- 
function perturbations. Let be the probability for any temperature pixel 
to be perturbed i units at the jth Hubble time-step where i = —Q, ...,+Q 
and j = 0, ...,M. We now have a total of 2Q + 1 possible outcomes for 
each 'trial'. Therefore, the distribution function can be obtained from a 
generalization of the trinomial: the multinomial distribution. The multino- 
mial distribution fn{x{l), ...,x{R)) gives the probability that an experiment 
consisting of n 'trials' each with R possible outcomes will result to x{i) 
occurrences of the ith outcome (i = 1,...,R) given that the probabilities 
for each outcome are p{l), ...,p{R). In direct correspondence with the tri- 
nomial, fn{x{l), ...,x{R)) is the general term of the multinomial expansion 
(p(l) + p{2) + ... +p{R))^. The interesting variable for our purposes is the 
total temperature fluctuation which at the jth time-step is given by 



kj = ^ i(x;- — X 



i=l 



where x'j is the number of i unit perturbations at the jth Hubble time-step. 
Using the multinomial expansion and equation (3.3), it is straightforward to 
obtain the characteristic function for the distribution of kj. The result is: 

Q Q 
cj){nj,p],...,pf,uj) = {2Y,p]cos{i uj) + (1 -2Y,p]r^ (3.12) 

1=1 i=l 

In order to proceed further we must specify the probability distribution of 
the step-function magnitudes i.e. the dependence of p*- on the index i. The 
simplest physically interesting choice is the distribution 

2 X 4^^-i 

which corresponds to a uniform distribution of the parameter /? in equation 
(2.1) in the range [0, Pmax] where /dmax should be chosen to be about unity. 



= 7^ (3-13) 
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Using now (3.12), (3.13) and (3.10) in (3.8) we obtain the generalization 
of (3.11) appropriate for cosmic strings: 



M~16 Aj-M Q 

^s{^) = n i-TT E + (1 - (3.14) 

j=0 ^ i=l 

which along with (3.11) is the central result of this work. In what follows 
we use the results (3.11) and (3.14) to obtain the probability distributions 
and several moments for temperature fluctuations induced by textures and 
cosmic strings. In accordance with the accepted values of the scaling so- 
lution parameters we will use rih = 10 for strings|p|[^ and rih = 0.04 for 



textures |30[]. Our goal is to compare the derived results with those cor- 
responding to a Gaussian distribution. It is therefore instructive to first 
investigate if and under what conditions can we obtain a Gaussian limit for 
the characteristic functions (3.11) and (3.14). 

Consider first the texture case described by equation (3.11). Since we 
want to compare with the standard Gaussian distribution which has a = 1 
we must first appropriately normalize the variable K dividing by its standard 
deviation to match the standard deviation u = 1 of the standard Gaussian. 
Consider the new variable 

K„ = y—^L== (3.15) 
The characteristic function that generates the moments of Kg is 

y/nh{M + 1) 
It is now straightforward to show that 

j=0 

But e~'^^/^ is the characteristic function for the standard Gaussian p9[| distri- 
bution. Therefore the distribution of the appropriately normalized variable 
Kg approaches the standard Gaussian for rih — > oo. Notice that the Gaus- 
sian limit is not obtained for M — > oo i.e. for the Gaussian limit to be 
realized we need several perturbations per horizon volume but evolution 
over more Hubble time-steps does not help. In a similar way, it may be 
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shown that the Gaussian hmit for the string characteristic function (equa- 
tion (3.14)) is obtained for the variable 



M ^ 

" § Vnh{M + l){Q + l){2Q + l)/6 ^^'^^^ 

provided that — > oo. For this variable the string characteristic function 
is: 

= ) (3.19) 

4 Results-Discussion 

We are now in position to use (3.11) and (3.14) (or (3.16) and (3.19)) to make 
predictions and compare these with the predictions of the corresponding 
Gaussian distribution. We first obtain the moments of the distributions. 
This may be achieved by differentiating the characteristic functions and 
using equation (3.5) to obtain the moments < K"^ >. Alternatively we may 
expand the characteristic functions around uj = 0. It is the later method we 



used here. We used the symbolic manipulation package Mathematica[31\ to 
expand (3.16) and (3.19) around a; = up to order 12. Prom the values of 
the coefficients and equation (3.5) we obtained the 12 lowest moments for 
the texture and string cases and also for the standardized Gaussian. The 
values of these moments are shown in Table 1. 



Table l:The values of the lowest six non- vanishing moments for the 
string (< >s), texture (< >t) and standard Gaussian (< Kjp >g)- 



m 


< K >s 


< >t 


< >G 


2 


1 


1 


1 


4 


3.0045 


4.124 


3 


6 


15.0675 


35.5577 


15 


8 


105.947 


440.226 


105 


10 


959.239 


9681.3 


945 


12 


10630 


159955 


10395 



The quantity of interest is the relative deviation from the Gaussian de- 
fined as: 



R„ 



< >G 



(4.1) 
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where < K^p >s^t is the standardized string or texture mth moment while 
< K"^ >G is the corresponding standard Gaussian moment. 

In Fig. 4 we plot the relative deviation Rf^ versus m for strings (notice 
that only even moments are shown since all odd moments vanish) obtained 
by expanding (3.19) with M = 16, Q = 5 and n/j = 10. Clearly, Rf^ is a 
rapidly increasing function of m which is evidence for the presence of long 
non-Gaussian tails in the string distribution function. However, even for the 
m = 12 moment the relative deviation does not exceed 3% implying that 
the non-Gaussian features in the string distribution function are fairly weak. 
The reason for this is that in the case of strings we have n/j = 10 ^ 1 which 
implies that the Gaussian limit is approached in an effective way. We have 
performed tests for different values of Q and found that the values of i?^ 
remain insensitive to within a factor of 2. 

In Fig. 5 we plot the corresponding relative deviations i?^ for textures. 
We used equation (3.16) with M = 16 and = 0.04. In the case of 
textures, not only is R^^ a rapidly increasing function of m but also even 
the lower moments are significantly larger than the corresponding Gaussian. 
For example the kurtosis (defined as < Kg > /(< Kg > )) is predicted to 
be 40% larger than the kurtosis of the standard Gaussian distribution, while 
the sixth moment is larger by a factor of three. As in the case of strings, the 
skewness and all the odd moments are found to vanish. This is due to the fact 
that the superimposed kernels are symmetric with respect to positive and 
negative perturbations. Clearly, such an assumption even though reasonable 
for CMB consideration is inapplicable for large scale structure calculations 
where a non-zero skewness is predicted by seed-based models [p^] [pO|]. 

The characteristic functions (3.11) and (3.14) can also be used to find 
the temperature fiuctuation distribution functions. These can be obtained 
by Fourier transforming the characteristic functions according to equation 
(3.4). The results may then be compared with the corresponding Gaussian 
with the same standard deviation (in order to keep the Fourier transform 
simple, we use the original forms (3.11) and (3.14) in this case). In Fig. 
6a we show the distribution function Fs{K) for strings, obtained by Fourier 
transforming (3.14) with M = 16, Q = 5 and n/j = 10. The difference 
Fs{K) — Fg{K) between the string induced distribution function Fs{K) 
and the corresponding Gaussian Fg{K) is shown in Fig. 6b. The relative 
difference {Fs{K) — Fg{K)) / Fg{K) at any given point does not exceed 1% 
but the presence of long non-Gaussian tails is clear. It is these tails that 
cause the rapid increase of the moments with the order m. 

Fig. 7a shows the distribution function for textures obtained by Fourier 
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transforming (3.11) with M = 16 and Uh = 0.04. Superimposed is the 
corresponding Gaussian distribution function. In this expected 
since n/^ <C 1, the non-Gaussian features are fairly clear. The central 
pronounced peak and the long tails seem to be generic features for tem- 
perature fluctuations induced by topological defects. Fig. 7b shows the 
difference Ft{K) — Fg{K) between the texture distribution and the corre- 
sponding Gaussian. It shows the same features as Fig. 6b which applies to 
strings but in the case of textures the magnitude of the relative difference 
{Ft{K) — Fg{K))/ Fg{K) is almost two orders of magnitude larger. This is 
reflected in the relative deviation of moments for which we have ^ i?^. 

Our results on both the relative deviation of moments from the Gaussian 
and the distribution function itself show the following: 

1. The non-Gaussian signature of cosmic strings is difficult to detect by 
measuring relative deviation of moments from the Gaussian. Relative 
deviations need to be measured to within less than 1% in order to 
distinguish cosmic strings fluctuations from Gaussian. The origin of 
this approximately Gaussian behavior of string induced perturbations 
is the large number of strings per horizon volume {n^ ~ 10). The 
large number of superimposed non-Gaussian features 'averages out' to 
an approximately Gaussian result as predicted by the central limit 
theorem. 

2. The measurment of moments provides a much more powerfull test for 
the texture model. This is due to the small number of textures un- 
winding per horizon volume {rih ~ 0.04) which avoids the suppression 
of the texture non-Gaussian features. Measuring the relative deviation 
of the kurtosis to within 40% should be enough to detect the deviation 
induced by textures. For the relative deviation of the sixth moment, 
a measurment accurate to within a factor of less than three, would be 
sufficient to indicate the presence of a texture signature. 

The tests based on relative deviation of moments from the Gaussian that 
have been studied here, have several interesting and powerful features, par- 
ticularly for testing the texture model. However, they are not sensitive to 
geometrical and topological features of the temperature fluctuation maps. 
Such features have been examined in Ref. |21] using a numerically obtained 



realization of cosmic string induced perturbations. It was found that topo- 
logical and geometrical tests can be a sensitive probe of stringy non-Gaussian 
features. An interesting extension of the work presented here is the study of 
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the geometrical features of string and texture induced temperature patterns 
using analytical methods and Monte Carlo simulations. Such a project is 
currently in progress j32[|. 
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6 Figure Captions 

Figure 1: The production of step-like discontinuities in the microwave 
temperature for photons passing on different sides of a cosmic string S. The 
string deficit angle is a and O is the observer. 

Figure 2: The propagation of a photon beam from the recombination 
time tree to the present time to- The horizon in three successive Hubble 
timesteps is also shown. 

Figure 3: The effects of a step-function perturbation on an initially 
uniform one dimensional distribution. 

Figure 4: The relative deviation of moments from the standard Gaus- 
sian. Rf^ corresponds to moments due to string induced perturbations and 
is plotted versus the m where m is the order of the moments. Odd moments 
are omitted since they vanish. 

Figure 5: The relative deviation i?^ for the case of textures. The de- 
viations from the Gaussian are much larger compared to the case of strings. 

Figure 6a: The distribution function Fs{K) for string induced pertur- 
bations. 

Figure 6b: The difference Fs{k) — Fg{K) where Fg{K) is the Gaussian 
distribution with the same standard deviation as Fs{K). The relative dif- 
ference does not exceed 1% but the presence of long non-Gaussian tails is 
clear. 
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Figure 7a: The distribution function Ft (K) for texture induced pertur- 
bations superimposed with the Gaussian distribution of the same standard 
deviation. 

Figure 7b: The difference Ft{k) — Fg{K) where Fg{K) is the Gaussian 
distribution with the same standard deviation as Ff [K] . The relative differ- 
ence exceeds 10% and is much more prominent than in the case of strings. 
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